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The propagation of reconnection signatures and their associated energy are examined using kinetic 
particle-in-cell simulations and Cluster satellite observations. It is found that the quadrupolar out-of- 
plane magnetic field near the separatrices is associated with a kinetic Alfven wave. For magnetotail 
parameters, the parallel propagation of this wave is super-Alfvenic (Vji ~ 1500 — 5500 km/s) and 
generates substantial Poynting flux (S ~ 10~ 5 — 10~ 4 W/m 2 ) consistent with Cluster observations 
of magnetic reconnection. This Poynting flux substantially exceeds that due to frozen-in ion bulk 
outflows and is sufficient to generate white light aurora in the Earth's ionosphere. 

PACS numbers: Valid PACS appear here 



Magnetic reconnection plays an important role in many 
plasma systems by releasing large amounts of magnetic 
energy through the breaking and reforming of magnetic 
field lines (e.g., [l|). During magnetospheric substorms 
the global magnetic geometry is reconfigured @, releasing 
magnetotail magnetic energy and creating intense auro- 
ral disturbances. During solar flares, magnetic energy 
release in the corona energizes large numbers of electrons 
which create hard x-rays when they impact to surface of 
the sun (e.g., ||). 

The sudden onset of magnetospheric substorms is be- 
lieved to be caused by either a near Earth instability at 
around 10i? e downtail (e.g^, [H) or reconnection onset 
around 20 to 30i? e (e.g., |5j|). Determining the mech- 
anism or mechanisms which are most relevant requires 
careful timing studies and has been the subject of much 
scrutiny and controversy (e.g., [3-Uoj). A key unanswered 
question regarding magnetic reconnection, therefore, re- 
gards how fast the released energy and associated sig- 
natures propagate away from the X-line. The propaga- 
tion of MHD signatures, ion flows and magnetic distur- 
bances, has been extensively studied in both substorms 
(e.g., [TTj] ) and solar flares (e.g., [13), but these mecha- 
nisms are limited by the Alfven speed. In some substorm 
events, however, it has been reported that the time lag 
between reconnection onset and auroral onset was less 
than the Alfven transit time from the reconnection site 
to the ionosphere (rj [l3j. It is necessary, therefore, to de- 
termine the nature of the reconnection signal that prop- 
agates fastest away from a reconnection site and its asso- 
ciated energies. Poynting flux[l4|, UH associated with ki- 
netic Alfven waves, for example, has been postulated as a 
possible energy source for aurora fi~6f . with observations of 
these waves near magnetotail reconnection sites (l7l [l8j. 

We simulate magnetic reconnection with the kinetic 
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particle-in-cell code P3D and find that the quadrupolar 
Hall out-of-planc magnetic field located near the separa- 
trices is associated with a kinetic Alfven wave (KAW). 
This KAW magnetic field perturbation has a supcr- 
Alfvenic parallel propagation speed (using lobe densi- 
ties), and is associated with a substantial Poynting flux 
that points away from the X-line. This KAW will ex- 
ist whenever Hall physics is active in the diffusion region 
fl9j . Simulation Poynting flux is consistent with Cluster 
statistical observations of multiple magnetotail reconnec- 
tion events. Scaling to magnetotail and ionospheric pa- 
rameters, the transit time of this standing KAW from a 
near Earth X-line is on the order of 50 seconds. 

Simulations: Our simulations are performed with the 
particle-in-cell code p3d t 20|. The results are presented 
in normalized units: the magnetic field to the asymp- 
totic value of the reversed field Bo, the density (no) to the 
value at the center of the current sheet minus the uniform 
background density, velocities to the Alfven speed ca, 
lengths to the ion incrtial length di, times to the inverse 
ion cyclotron frequency fi" 1 , temperatures to m^c^, and 
Poynting flux to So = caBo/^tt. We consider a system 
periodic in the x — z plane where flow into and away 
from the X-line are parallel to z and x, respectively. The 
initial equilibrium consists of two Harris current sheets 
superimposed on an ambient population with a uniform 
density of 0.2. The equilibrium magnetic field is given by 
B x = tanh[(i/ - L z /4)/w ] - tanh[(y - 3L z /4)/w } - 1, 
where wo and L z are the half- width of the initial current 
sheets and the box size. The electron and ion temper- 
atures, T e = 1/12 and Tj = 5/12, are initially uniform. 
The simulations presented here are two-dimensional, i.e., 
d/dy = 0. Reconnection is initiated with a small initial 
magnetic perturbation. 

We have explored the separatrix structure and recon- 
nection signal with three different simulations, varying 
the electron mass as shown in Table [TJ These simula- 
tions were used for a previous study of reconnection pTj. 
As shown in Fig. 1 of citation[2l| the reconnection rate 
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TABLE I: Simulation parameters and results: A = grid scale, 
c = light speed, (L X ,L Z ) = system size, (A, k, V sim , S sim ) = 
properties of KAW. (du,d e e) = lobe inertial lengths. 
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FIG. 1: 2D overview plots, (a) B y , with red and blue mag- 
netic field line segments used in Fig. [2] (b) Parallel poynting 
flux S • h xz with ion velocity vectors. Black line is separatrix 
magnetic field line. 



increases with time, sometimes undergoes a modest over- 
shoot, and approaches a quasi-steady rate of around 0.15. 

The structure of roughly one quadrant of the recon- 
nection region is shown in Fig. [TJ with the X-line lo- 
cated at (x/d i7 z/di) = (21.36,6.40). The large B y as- 
sociated with Hall physics is clearly evident near the 
separatrix. This magnetic field is produced by nearly 
parallel electron flows near the separatrices, which are 
strongly super- Alfvenic. There is a strong Poynting flux 
parallel to the in-plane magnetic field, S • b X2 , where 
b xz = (B x + & z )/yjBl + B 2 . It is this Poynting flux 
which carries the energy of the first signal of reconnec- 
tion. Note that there is little ion flow associated with 
this By and Poynting flux. 

In order to gain some handle on the physics governing 
this By structure associated with reconnection, we rep- 
resent it as a superposition of linear waves with various 
k values. The scaling laws based on this analysis will be 
shown to be consistent with simulation properties. Ex- 
amining Fig.[T]a,, the quasi-lD B y structure is very nearly 
parallel to the separatrix and thus is a strongly oblique 
wave with fc|| <C k±. As a starting point, we use the two- 
fluid analysis from previous studies pjl, [22|, and analyze 
the branch of waves associated with Alfven waves and ki- 
netic Alfven waves. For the simulation parameters used 
in this study and noting also that kd e < 1, with d e the 
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with D = 1 + k 2 d\ and c 2 = (T e + T^/m,. Note that for 
highly oblique waves, the parallel group velocity is equal 
to the parallel phase velocity 

Simulation KAW: This analysis uses quasi-steady 
reconnection to study the properties of the separatrix 
kinetic Alfven wave (KAW). This is necessary because 
the KAW is so fast that the quadrupolar field associated 
with it very quickly fills the whole simulation domain 
making velocity measurements due to direct time vari- 
ation impossible; the high speed of the KAW is most 
likely why its propagation velocities have been largely 
ignored by previous studies of collisionless reconnection, 
although other properties of the quadrupolar field have 
been extensively examined through simulations, satellite 
observations, and laboratory experiments ([l[, and ref- 
erences therein). During steady reconnection, magnetic 
field lines convect along the inflow (z) direction, recon- 
nect, and then flow outwards. The propagation veloc- 
ity of the KAW can be measured by changing frames 
to one moving with that inflowing magnetic field line. 
Since the reconnection is steady, the time difference be- 
tween two magnetic field lines is the difference in flux 
between the two lines over the reconnection rate, i.e., 
At = Atp/E r = Ai/j/(d/dt(ip xU nc - Valine)), where ip is 
defined such that B = y x Wtp + B y y. By examining 
the KAW By at different ip values, therefore, one can de- 
termine the propagation speed of the KAW. An example 
of this analysis for the m e /irii = 1/400 case is shown in 
Fig. [2 which shows the variation of B y along magnetic 
field lines (lines of constant ip). For clarity, two repre- 
sentative magnetic field line segments colored red and 
blue arc shown in Fig. QJi, and the B y plots taken along 
them are colored the same. Each B y plot represents a 
Aip = 0.05 Bodi, and each successive plot has been off- 
set 0.08 Bq along the vertical from the previous one. The 
evolution of the B y is not characterized by a simple prop- 
agation. First, the peak value of B y increases with time. 
Second, the dispersive nature of KAWs also leads to mul- 
tiple velocities associated with the B y structure. The lo- 
cation where B y = propagates at the peak KAW speed, 
Vpcak ~ C se = (T e + Ti)/m e w 14 ca- However, there 
is little Poynting flux associated with this velocity. In- 
stead, we focus on the propagation of the main B y signal 
by finding the velocity of the wave front. The two dashed 
lines in Fig. [2] have the same slope and denote the wave 
front in two of the curves separated by Aip = 0.15 £>o<i- 
The propagation velocity of the x-intercept of this slope 
(shown as vertical green lines) is calculated to be 4.0 ca, 
which is substantially less than the peak parallel KAW 
speed. The measured values are shown as V S i m /cA in 
Table H 

It is critical to determine if this propagation velocity 
is consistent with the kinetic Alfven wave predictions of 
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FIG. 2: By plotted along magnetic field lines separated by 
Alp — 0.05 Bodi, which represents At = 0.31 . Vertical 
offset of each plot is 0.08 -Bo- Wave fronts for Aif) = 0.2 Bodi 
and 0.35 Bo^i shown as dashed lines, with respective B y — 
shown as horizontal dotted lines. Blue and red B y plots 
taken along field line segments shown in Fig. Wave front 
intersections with B y = denoted with vertical green lines. 



Eq. (Q}. First, the k values associated with this B y Hall 
field must be determined. Vertical slices of the Poynting 
flux S ■ h xz were analyzed at the locations of the wave 
front (m 4 /m e = [25,100,400], x/d t = [170.0,90.0,35.0]). 
The magnitude of this Poynting flux is shown in Table U 
as S s i m . The width at half- max of the Poynting flux was 
measured and used to determine the primary k = 2n/X 
value for the KAW. As an example, for the mi/m e = 
400 case, the half max was 5 = 0.65 d,;, yielding A « 
2.6 di. The standing KAW wave is located close to the 
scparatrices, so the simulation lobe plasma values are 
used to determine parameters (B w 1.0, n «a 0.2, Ti + 
T e 0.5), which yields kdu 5.4, where the denotes 
lobe values. The resulting \,kdu, and kd e £ are shown 
in Table |U Plotting the velocities predicted from Eq. ^ 
versus the simulation measured velocities yields excellent 
agreement, as shown in Fig. 

Associated with this Hall structure are electron beams 
and significant Poynting flux. The super-Alfvenic elec- 
tron beams are associated with the parallel currents 
which create the quadrupolar B y . A theoretical predic- 
tion for the Poynting flux can be determined for com- 
parison with simulation values. We use S ■ h xz w S x = 
(c/4tt)(E x B) x w -(c/An)E z B y . The normal Hall elec- 
tric field is due to the frozen-in electron flow, which 
dominates over the ion flow, giving E z ~ Vex-E$y f c ~ 
-J x B y /(nec) « -B y B' y /(Aime) with B' y = dB y /dz. 
Substituting gives S x sa B y B' CAy di/Air, where ca v = 
Byj \J A-nrriin. Note that the integrated KAW Poynting 
flux is independent of the width of the KAW. As with the 
KAW velocity determination, n is the lobe density with 
B y si 0.25 consistent with simulation values. Comparison 
of the theoretical Poynting fluxes with simulation values 
also yields excellent agreement, as seen in Fig. [SJd. Note 
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FIG. 3: Simulation values versus theoretical predictions for 
(a) KAW propagation speed and (b) Poynting flux. 



that this KAW Poynting flux substantially exceeds the 
Poynting flux associated with the ion bulk flow away from 
the X-line: S x , ion c a B 2 /Att, since B 2 Z 0.01 B 2 <C B*. 

Comparisons with Satellite Data: A statisti- 
cal study of reconnection events has been performed 
previously [2^], where magnetotail reconnection cross- 
ings with correlated Geocentric Solar Magnetospheric 
(GSM) B z and Vi X reversals were selected. In that 
study [23[, comparisons with simulations were made by 
renormalizing data using magnetic fields just upstream 
of the separatrices (B s ) and densities in the ion out- 
flow region (ri ou t), yielding normalization velocity ca = 
B s /-iJ 47rm,Tiout an d Poynting flux S = Ga B 2 /Air. Using 
these normalizations, the Poynting flux from this Clus- 
ter data set is compared with data from the mi/m e = 25 
case. For the simulation data, the normalization values 
used were B s = 0.8 and n out = 0.2. The simulation 
sub-region used was a rectangle roughly centered on the 
X-line with length approximately 35 di and height ap- 
proximately 13 di, using n = 0.2. Fig. U shows this com- 
parison, where only normalized S x /S > 0.02 is plotted. 
Tailward S x /S is shown in red and Earthward S x /S is 
shown in black. The bounds of the simulation and Clus- 
ter data are similar, being limited to \Vi X \ ~ 0.7 and 
\B X \ ;$ 1.0. The separatrix KAW structure is present in 
both plots in the region of large B x and nearly zero Vi X . 
Both datasets show a strong correlation in the sign of Vi X 
and S x , implying that the Poynting flux points away from 
the X-line. However, for small Vi X and larger B x there is 
some anti-correlation which corresponds to ion flow to- 
wards the X-line just outside the separatrices. Both data 
sets show significant S x for small \B X \ and larger negative 
Vi X , which is associated with the very long outflow jet of 
super-Alfvenic electrons seen in simulations with kinetic 
electrons [U and satellite observations [25| . There is 
an asymmetry, however, in the satellite data along Vi X 
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FIG. 4: (a) Simulation and (b) Cluster Observations: Poynt- 
ing flux S x /S > 0.02 scatter plot in (B x /B a ,Vi x /cA) plane 
(red: Tailward S x /S, black: Earthward S x /S). Simulation 
values in (a) are smaller circles due to large number of data 
points. Normalizations described in text. 



not present in the simulations, with only negative (tail- 
ward) S x /S having significant values for B x « and 
finite Vi X . Some possible explanations are: (1) In most 
of the events, the satellite was initially tailward of the 
X-line and then crossed to the Earthward side, so Earth- 
ward flows represent more developed X-lines. (2) The ob- 
stacle presented by the strong Earth's dipole field could 
create back pressure and lead to outflow asymmetries at 



the X-line. Or (3) 3D effects lead to this asymmetry. 

Predictions for the Magnetotail: The KAW as- 
sociated with the quadrupolar B y propagates at a super- 
Alfvcnic speed and carries significant Poynting flux. To 
assess its importance for the magnctosphcre, we use 
the following typical parametersjg: B w 20 nT, n « 
0.1cm- 3 , T e w 300 eV, and T, IkeV. As the KAW 
propagates large distances in the magnetotail, it is quite 
probable that the k associated with it will decrease ow- 
ing to the dispersive nature of KAWs. Taking the sim- 
ulation kd e « 0.3 to be the maximum expected k, we 
take kdi w 1 to be the minimum k because at this k the 
KAWs are no longer dispersive. As is found in the simu- 
lations, we use By/Bi \, c « 0.25. These values yield the 
following ranges of parameters associated with the KAW: 
Vj| ~ 1500 - 5500 km/s, S ~ 0.7 • 10~ 5 - 9 • 10~ 5 W/m 2 . 
For an X-line located 20 R e downtail from the Earth, the 
predicted propagation time is At ~ 25 — 85 sec, which is 
substantially less than the Alfven transit time (~ 250 sec) 
for the same distance. 

An important question remains as to whether this 
KAW energy will be able to propagate to the Earth's 
ionosphere and create aurora. In the simulations (largest 
L x w 10 R e and At ss 50 s using simulation lobe parame- 
ters), the KAW propagates all the way to the edge of the 
simulation, but the limited length scale as well of lack 
of a dipole geometry make exact estimation of the wave 
attenuation impossible. This is an important question 
currently under study. Assuming parallel propagation of 
the Poynting flux so that it stays on the same magnetic 
flux tube, the Poynting flux in the ionosphere Sion would 
be: Sion — (Sion/Siobc)5'iobo ~ 10 3 Siobc- Reducing this 
flux by a factor of ten as an estimate of attenuation yields: 
Sion ~ 10 2 5iobc ~ 0.7 • 10~ 3 - 9 • 10- 3 W/m 2 , which is 
still on the order of or greater than the 10 _3 W/m 2 = 
1 ergs/cm 2 s necessary to create a white light aurora. 
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